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ABSTRACT 


In  this  paper  we  introduce  a  new  adap¬ 
tive  processor  able  to  work  well  even 
when  the  desired  signal  and  the 
interference  are  coherent.  The  present 
adaptive  processors  tail  to  operate  in 
these  cases.  The  results  of  simulations 
appear  to  confirm  the  theoretical 
predictions. 

I.  Introduction 

Since  the  pioneering  work  of  Howells[tj  , 
Applebaum[2]  and  Widrow[3j  ,  there  has  been  con¬ 
siderable  activity  in  the  development  of  adaptive 
beamforming  for  radar,  sonar,  communication, 
spectral  estimation,  efc  . 

Though  the  details  differ  in  the  different  appli¬ 
cations  the  main  assumptions  and  processing  algo¬ 
rithms  are  essentially  the  same.  In  particular,  a 
key  assumption  in  all  the  previously  cited  work  is 
that  the  interfering  signals  are  not  coherent  with 
(i.e.,  do  not  have  fixed  phase  differences  from)  the 
desired  signal.  More  generally,  two  signal  will  be 
said  to  be  coherent  if  one  is  a  scaled  and  delayed 
replica  of  the  other.  Coherent  interference  can 
arise  when  multipath  propagation  is  present,  or 
when  "smart"  jammers  deliberately  introduce 
coherent  interference,  e  g.  by  retrodirccting  the 
signal  energy  to  the  receiver.  Coherence  can  com¬ 
pletely  destroy  the  performance  of  adaptive  array 
systems. 

In  the  referenced,  5]  Gabriel  and  Evans  et  of 
show  that  the  subaperture  sampling  or  spatial 
smoothing  idea,  as  they  call  it,  can  be  applied  to 
the  off-line  eigenstructure  based  method  for  direc¬ 
tion  finding.  This  is  an  important  contribution,  the 
main  idea  of  which  we  independently  rediscovered 
later  (see  Shan  e<  al.  [6]  ).  However,  methods  of 
applying  this  idea  to  adaptive  beamforming  (rather 
than  direction  finding)  were  not  obvious.  In  this 
paper  we  will  introduce  an  adaptive  processor 
based  on  spatial  smoothing  algorithm. 
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II.  The  Optimum  Weight  Vector  of 
Adaptive  Processor  for 
Noncoherent  and  Coherent  Inputs 

It  is  easy  to  show  that  the  optimum  weight  vec¬ 
tor  for  different  criteria  is  of  the  form 

wopt  =  aR,"ja(T0) 

where  R„  =  Ex(t)x'(t)  is  the  covariance  matrix 
of  the  array  measurement  vector  x  ,  a  is  a  scalar 
constant,  and  a(r0)  is  so-called  "look"  direction 
vector. 

We  shall  assume  that  we  are  interested  in  the 
signal  s(<)  in  the  known  look  direction  t?D.  and 
that  the  interfering  signals  are  from  the  other  K-l 
unknown  directions  ji5, . tJa-i! 

Noncoherent  Signals 

Under  the  assumptions  of  statistical  indepen¬ 
dence  and  of  noncoherence,  we  shall  express  the 
array  measurement  data  vector  as 

x(f )  =  a(r0)s  (t)  +  7tj(f )  +  v(<) 

where  K  specifies  the  interference  directions 
X  =  (a(Ti)  ■  ■  •  a(r/f_,)] 

and  s(f)  ,  the  array  input  signal  vector,  is  consist¬ 
ing  of  desired  signal  s  (f )  and  jamming  signals  j  (f ) 


s(<) 

s  (t  ) 

S«)  = 

Ji(t ) 

j(0 

7at-i(0 

The  covariance  of  x(f)  can  now  be  written 
R„  =  pf  a(T0)a*(r0)  +  Rnn 


where 


Rim  =  +  o2\ 


and  we  have  assumed,  for  simplicity,  the  noise 
intensity,  a2,  to  be  same  at  each  sensor 

By  using  the  matrix  inversion  lemma,  we  can 
show  that 


w  opt 


=  0 
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where  0  is  a  scalar  constant 

We  also  introduce  the  modal  representation 

=  E  Keiei  . 


As  =  a(r0)s(f)  +  £  a(rjy,(<) 

1*1 

=  [a(r0)  +  7,a(T,)  +  +  7K-MrK-u'.s  (t ) 


where  S\j  and  { e*  {  are  the  nonzero  eigenvalues 
and  the  corresponding  eigenvectors  of  the  M  x  M 
matrix  XR^X*.  which  will  have  rank  K  -  1 
because  RJ;  is  the  covariance  matrix  of  the  K  —  1 
noncoherent  signals. 

Finally  we  shall  also  assume  that  the  back¬ 
ground  measurement  noise  intensity  is  small  com¬ 
pared  to  the  signals  y(f).  so  that  we  shall  have 

\  »  a2 


where  the  jy,  |  are  fixed  complex  constants  given 
by 

7t  =  (p,/Po>  ,  i  =  1 . K-\ 


In  this  case,  the  covariance  matrix,  AEsa’A',  will 
have  rank  1.  so  that  it  will  have  only  one  nonzero 
eigenvalue  A],  and  the  covariance  matrix  RZI  will 
have  M  -  1  eigenvalues  equal  to  o2. 

Therefore  we  shall  have 
"opt  =  aRzi' a(r0' 


\  +  (T2 


Then  we  can  write 
"opt  =  0Rnn-1“('r o) 


K- 1  ,  U  .  J 

E  - 5“ ® t  ®  t*  *  E  2~ ®  t  •  t*  »(To) 

t  =  i  Al  +  a  i=jr  o' 

o  tf 

=  2  S  P.e*  •  Pi  =  eja(T0)  . 


°z  i~K 


where  =  denotes  asymptotic  (  as  o2-*0  )  equality. 
Now  by  construction,  the  direction  vectors 

|a(r,) . a(T*_,)j  of  the  interfering  signals,  which 

are  the  columns  of  the  matrix  X,  lie  in  the  span  of 

the  first  K  —  l  eigenvectors  je, . e*.,!  and 

are  therefore  orthogonal  fo  the  remaining  eigen¬ 
vectors  [e  k . e  j|f ). 

Therefore,  we  shall  have 

*op*a(Ti)  =  Piei“(n)  -01  =  1 . K- 1  . 

cr  t=/r 


so  that  the  beam  pattern  will  have  "deep  nulls"  in 
the  interference  directions.  In  the  look  direction, 
on  the  other  hand,  the  constraint  will  ensure  that 
we  have 

"  opt  a  (  tq)  —  1 

However,  the  situation  deteriorates  badly  in 
the  coherent  case. 

Coherent  Signals 

If  the  signals  have  fixed  phase  differences, 
which  rally  means  equal  frequencies  and  fixed  i px, 
then  we  shall  have  the  representation 


f  j  u  , 

-a  — — —  e  i e  ’  +  V  — - e.e,*  a(r0) 
A,  +  o-  i-2  cr 


=  a  V  -^e/aCo^e, 

t=z  a- 

All  wo  can  aay  hare  ib  that  the  linear  combination 
b  :=  a(r0)  +  y2a(T,)  +  •  ■  •  +  y*-a(T*-_,) 

will  be  orthogonal  to  the  \e2 . eu\.  This  does 

not.  however,  imply  that  the  same  will  be  true  of 

the  ia(T)) . a(T*_,;j  individually,  and  therefore 

there  will  not  in  general  be  any  nulls  in  the  direc¬ 
tions  of  the  interfering  signals. 

Due  to  the  constraint,  we  shall  have 

"opt*(T0)  =  1  , 

but  this  is  of  small  comfort,  because  the  actual 
array  output  will  be 

y  (0  =  w„'p(x(f) 


I  U  ,  r 

a|E  -V(cta(T,))e,'  bs(0  +  v(0 


where  we  recall  that  b  :  =  a(T0)  +  V  7ia(T,)  lies 
along  ej  and  is  orthogonal  to  |ej . Cyj 

Therefore  there  will  be  no  signal  output  from  the 

conventional  array  when  the  signals  are  coherent 

A  Way  Out 

This  analysis  also  makes  clear  what  is  neces¬ 
sary  to  rescue  the  situation:  we  must  somehow 
restore  the  rank  of  the  covariance  matrix 
E{As)(As)’  to  being  K.  Then  the  noise-alone 
eigenvectors  will  be  orthogonal  to  all  the  vectors  in 
the  space  of  the  signals  and  the  beam  pattern  will 
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have  nulls  in  the  directions  of  the  interfering  sig¬ 
nals.  A  simple  scheme  for  achieving  this  rank  res¬ 
toration  with  an  adaptive  algorithm  is  proposed  in 
the  next  section. 


III.  A  New  Adaptive  Array  Processor 

We  shall  describe  a  preprocessing  scheme  for 
the  sensor  outputs  that  will  restore  the  rank  of  the 
input  signal  covariance  matrix  to  K  even  if  the 
inputs  are  completely  coherent  with  each  other. 

The  scheme  is  based  on  combining  measure¬ 
ments  from  overlapping  subarrays. 

Given  the  M  sensor  outputs  at  any  time 
instant. 

x(<)  =  [xjCO  •  ■  •  xu(t)]T  . 

define  p  subsets  (recall  that  K  is  the  number  of 
sources' 

z(1)(0  =  [x,(0  •  •  ■ 

zt2>(0  =  [x2(t)  •  •  ■  **+a(t):r 


z (p,(<)  =  [xp  ■  ■  XK^f(t)'T  . 

Define 

r£>  =  **<*>*»>• 

and  the  spatial  smoothed  correlation  matrix 

r  = 

Pk* 1 

Then  we  can  prove  (  see  Referencef?,  0]  '  that  R 
will  have  the  form 

R  =  ASA*  +  ozI 

where 

S  will  have  rank  K  if  and  only  if  p  a  K  . 

Once  S  has  rank  K,  then  the  noise  eigenvec¬ 
tors  will  be  orthogonal  to  the  columns  of  A  and  by 
the  analysis  of  Section  III.  will  give  nulls  in  the 
interference  directions. 

Fig  la  shows  how  to  form  the  spatial  data  sub¬ 
set  from  one  snapshot'  A  flow  diagram  of  the  pro¬ 
posed  adaptive  processor  is  shown  in  Fig  lb.  At 
each  time  instant,  the  'snapshot'  of  M  data  samples 
is  divided  into  overlapping  subgroups  of  K  samples 
each,  these  subgroups  arc  then  fed  in  succession 
into  the  adaptive  processor,  which  updater  a  K- 
dirncnsional  (if  each  sensor  is  followed  by  a  N 
lappcd-dclay-lme,  then  a  A'xjV  -dimensional', 
weight  vector  each  time.  After  all  the  suogroups 


have  been  processed,  the  same  procedure  is 
repeated  with  the  next  ’snapshot'. 

It  is  important  to  note  that  the  processor  can 
be  any  kind  of  adaptive  processor  using  any  adap¬ 
tive  algorithm  and  any  array  structure,  e.g,  the 
Howells-Applebaum  or  the  Frost  arrays  could  be 
used.  There  are  several  ways  to  form  the  array  out¬ 
put  and  they  have  deferent  effects  on  the  perfor¬ 
mance  of  the  processor.  Further  studies  on  this 
subject  are  being  conducted. 


TV.  The  Computer  Simulation  Results 

In  the  computer  simulation  example,  we  have  a 
signal  0.1sin0.47Tf  arriving  at  90°  and  two  coherent 
interfering  signals  arriving  at  50°  and  130°.  Fig  2a 
shows  the  beampattern  of  a  conventional  Frost 
array  with  six  sensors.  Our  new  scheme  uses  10 
sensors  with  subarrav  of  size  six.  adapted  with  the 
same  LMS  algorithm,  and  gives  the  beampattern 
shown  in  Fig. 2b.  The  input  signal  waveform  is  shown 
in  Fig. 3a,  with  the  output  of  the  conventional  array 
in  Fig. 3b  and  that  of  new  processor  in  Fig  3c  - 
notice  the  big  dilTcrencc  in  the  scales  of  these 
figures.  The  results  would  seem  to  speak  them¬ 
selves 

Several  other  simulations  have  shown  similar 
results. 

V.  Conclusions 

Conventional  adaptive  processors  perform  very 
poorly  in  coherent  receiving  environments.  If  the 
received  signal  is  coherent  with  one  interference, 
the  signal  will  be  canceled  out  on  the  output,  the 
processor  will  therefore  totally  fail.  The  suggested 
new  adaptive  processing  system  is  able  to  over¬ 
come  this  degradation  of  performance  in  coherent 
input  environments,  without  considerably  increas¬ 
ing  the  complexity  of  the  system  structure  or  the 
computational  burden.  The  new  array  structure 
can  be  applied  in  conjunction  with  any  of  the  adap¬ 
tive  algorithms  and  structures  of  current  adaptive 
arrays,  and  successfully  separates  the  coherent 
array  inputs,  as  shown  by  theoretical  analysis  and 
simulation  results. 
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Figure  la:  Spatial  Smoothing 
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Figure  lb:  flow  Diagram  of 
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Figure  2a:  Beampattern  of  the  Frost 
Array  With  Coherent  Inputs 
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Figure  3b:  Output  of  Frost  Array  With 
Coherent  Interference 
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